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The phase diagram of the three-dimensional Z2 gauge Higgs system 
at zero and finite temperature* 

Luigi Genovese, Ferdinando Gliozzi, Antonio Rago and Christian Torrcro ^ 

'^Dipartimento di Fisica Teorica, Univcrsita di Torino and INFN, sezione di Torino, via P. Giuria, 1, 
1-10125 Torino, Italy. 

We study the effect of adding a matter field to the Z2 gauge model in three dimensions at zero and finite 
temperature. Up to a given value of the parameter regulating the coupling, the matter field produces a slight 
shift of the transition line without changing the universality class of the pure gauge theory, as seen by finite size 
scaling analysis or by comparison, in the finite temperature case, to exact formulas of conformal field theory. At 
zero temperature the critical line turns into a first-order transition. The fate of this kind of transition in the finite 
temperature case is discussed. 
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1. INTRODUCTION 

The three-dimensional Z2 gauge Higgs system 
is perhaps the simplest example of a gauge the- 
ory coupled to a matter field. Its action can be 
written as 

S = -Pg ^Ua- (3i^ (JxUxyOy (1) 

□ eA (xy) 

where [/□ = Ilfen^^' ~ ^'^^ A is a cu- 
bic lattice and cr = ±1 denotes the matter field. 
This model is self-dual under a Kramers- Wannier 
transformation: 

Z(/3g, Pi) = ^ e-^ cx Z{(3i. Pg) , (2) 

con/. 



with (3 — — log v^tajih/?. Its phase structure has 
been determined long ago ||l| and it has been 
shown to be very similar to that of SU (2) gauge 
system coupled to a matter field in the funda- 
mental representation, but of course it is much 
simpler, moreover the coupling to the Ising mat- 
ter can be now efficiently implemented by a non- 
local cluster algorithm which leads to very ac- 
curate Monte Carlo (MC) simulations. Therefore 
it appears as an ideal laboratory to test new ideas 
on the confining-deconfining properties of coupled 
gauge models. Recently it has been used to study 



the phenomenon of string breaking in order 

to probe a general mechanism proposed to explain 
why this phenomenon is not visible in the Wilson 
loops 1^. 

In this contribution we report an accurate anal- 
ysis of the transition lines of this model at zero 
and at finite temperature. 

At zero temperature we use the histogram 
method and the re- weighting techniques |^,^ 10 
to locate the first order transition line and the 
standard finite size scaling study of the continu- 
ous transitions, which turn out to be in the uni- 
versality class of 3D Ising model. We find that the 
apparent triple point suggested by the old anal- 
ysis based on the hysteresis cycles is a finite 
size effect. 

At finite temperature we argue that the matter 
field term in the action, for not too large values 
of is an irrelevant operator of the deconfin- 
ing transition of the pure Z2 gauge theory, which 
is known to belong to the 2D Ising universality 
class. We support this conjecture by comparing 
the Polyakov correlator with the exact formula of 
spin-spin correlator in a finite box given by the 
underlying 2D conformal theory. This allows us 
to locate the transition lines in the plane Pg,Pi ■ 
For large enough (]j the nature of the transitions 
changes and becomes strongly influenced by the 
boundary conditions. 



*Based on a poster by A. Rago and a talk by F. Gliozzi 
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Figure 1. The double peak structure of the his- 
togram of the hnk and the plaquette distribution 
near a first-order transition. 



2. ZERO TEMPERATURE 

The problem of detecting a first order transi- 
tion by MC simulations on a finite system of vol- 
ume V = L'^ can be solved by computing the his- 
togram of energy distribution P{E, L) at a point 
(3giPi close to the transition line and then ex- 
trapolating the data to nearby values |^,| 10 



P{E,L)^n{E,L) 



exp -PgE 
Z 



(3) 



where Z is the partition function and Q{E,L) is 
the number of states of energy E ^ — l^n ^ 
% J2{xy) o-xUxyO-y In the vicinity of a first or- 
der transition it has a characteristic double peak 
structure as shown in the double histogram of 
the links and plaquettes of Fig.(|l|). A suitable 
re- weighting through Eq.(||) yields the line /?/ = 
f{PG,L) where the two peaks at Ei{pG, Pi, L) 
and E2{Pg, Pi, L) are of equal height. We located 
numerically this line for cubic lattices of sides 
ranging from L = 10 to L = 30. The autocorrela- 
tion times for larger lattices were too large in our 
canonical MC simulations. Perhaps more refined 



techniques such as multi-canonical algorithm 
could reduce this correlation significantly. A typ- 
ical plot of Pg = f{Pi,L) is reported in Fig.(||): 
it is formed by two dual lines which cross each 
other on a self-dual point. Above the crossing 
point the distance of these two curves from the 
self-dual line (SDL) decreases rapidly when L in- 
creases and is already microscopic at L = 18 (see 
the inset of Fig.(H)). If the system had a true 
triple point, one should see a triple peak near the 
crossing point, while we observe in all the cases 
a sharp double peak structure as represented in 
Fig.(|l|). Denoting by {x,y) and {x',y') the coor- 
dinates of the two peaks in the (plaquette^ link) 
plane (see Fig.(||)) it is easy to prove that, if two 
states coexist at a self-dual point, they are re- 
lated, in the thermodynamic limit, by 



y' smh2Pi = cosh 2/3/ 
- y sinh 2/3/ — cosh 2/3/ 



(4) 

(5) 



which turn out to be approximately verified also 
in the finite lattices. It is worth observing that 
the presence of a double peak is not sufhcient to 
assure a true first order transition. A useful quan- 
tity in this regard is the bulk free-energy barrier 
AF between the two coexisting states, defined by 



AF{L) = W{E,r,)-W{E,) , 



(6) 
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Figure 3. Contour plot of the double peak struc- 
ture of Fig.(|l|) in the plaquette (x) and link (y) 
plane. 



where W = ~\ogP{E,L) and Em{{l3G,fii,L) is 
the local maximum which separates the two dips 
at El and E2 when W{Ei) = W{E2), as shown 
in Fig.(Q). At a continuous transition, AF(L) 
is independent of L and at a first-order transi- 
tion it increases monotonically with L. For large 
enough L one has |^ , for periodic boundary con- 
ditions, AF(L) ~ 2aL^, where a is the interface 
tension. In our data at fixed L AF is maximal 
at the crossing point. Extrapolating to large L 
we locate the point where the first-order transi- 
tion has its maximal strength at Pa = 0.708(2), 
Pi = — log(tanh(/3G))/2 and the corresponding 
interface tension is 2a — .00108(9). 

Using the behavior of AF{L) as a criterion for 
discriminating the order of the transition, we can 
prove the first-order nature only for a small in- 
terval around the crossing point. Near the bi- 
furcation, where the transition lines go off the 
SDL, even if small lattices show still the double 
peak structure, the transition is second order. To 
extract an estimate for the infinite-volume tran- 
sition line in this region, we tried the standard 
finite-size scaling form 

(3GiL) = f{l3j,^)-c{Pj)L-i, (7) 




1.6 10'' 1.7 Id* 1.8 Id* 1.9 irf 2.0 10* 2.110' 
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Figure 4. A typical plot of Vt^ = -logP(£',i) 
with i = 18 on the self-dual line. 



which turns out to fit well the data, as Fig.(||) 
shows, using v = 0.63, which is the value of the 
corresponding critical index of the 3D Ising uni- 
versality class. The resulting phase diagram in 
the thermodynamic limit is reported in Fig.(^. 

3. FINITE TEMPERATURE 

The universality class of a continuous deconfin- 
ing transition of a pure gauge theory in D dimen- 
sions is well understood in terms of the Svetitsky 
and Yaffe (SY) conjecture it coincides with 
that of the spin model in I? — 1 dimensions with a 
global symmetry coinciding with the center C{G) 
of the gauge group. What is the effect of adding 
matter to a pure gauge system? there is no gen- 
eral answer. Even when the matter can be treated 
as a 'small' perturbation ( using for instance the 
inverse mass 1/m as a perturbing parameter), it 
does not change the nature of the transition only 
if it is an irrelevant operator. 

At first sight one is tempted to conclude that 
the matter acts always as a relevant operator, 
since it breaks explicitly the global center sym- 
metry of the pure gauge theory at finite temper- 
ature, which is the hart of the SY conjecture. 

In the model at hand, when the coupling Pi of 
the matter is small enough, it is easy to show that 
this is not actually the case: the addition of the 
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Figure 5. Finite size scaling behavior of the tran- 
sition hne near the bifurcation. The three curves 
are the fits to Eq.(^. At each value of L the 
points are obtained by re- weighting the data of a 
single simulation. 



matter does not modify the universality class of 
the deconfining transition both at zero and at fi- 
nite temperature. The argument goes as follows. 
Performing a duality transformation as defined 
in Eq.(|^), when /3/ — > oo we recover the usual 3D 
Ising model. For large /3/ it is possible to do a per- 
turbation expansion in e^^' . The first order cor- 
rection is due to a single anti-ferromagnetic link, 
and the corresponding change in the free energy 



is proportional to e 



Near 



the critical point, this may be expanded in a sum 
of scaling operators all of which will be even un- 
der spin reversal. The dominant term is therefore 
proportional to the energy operator of the unper- 
turbed Ising model, both at zero and at finite T. 
Therefore the only effect is that of a slight shift of 
the transition line, without changing the univer- 
sal critical properties. So, in a sense, the matter 
field acts as an irrelevant perturbation of the uni- 
versality class of the pure gauge theory. In order 
to extend this property to larger values of /?/ we 
have to resort to numerical work. 

At finite T, where the deconfining transition is 
known to be well described by the 2D Ising uni- 
versality class, we can support this property even 
at larger values of /3/ by accurate numerical tests 



Figure 6. A schematic view of the phase diagram 
of the Z2 gauge Higgs model at zero temperature 
in the infinite volume limit. The cross indicates 
the point where the first-order transition has its 
maximal strength. 



of comparison with the exact finite size formulas 
dictated by the 2D conformal field theory. 

At /3/ = 0, A^t = 6 = 1/Tc the deconfining tran- 
sition is estimated to be at /3g = 0.746035 pl[]. 
Here the Polyakov loop correlators, according to 
the SY conjecture combined with the universal fi- 
nite size effects dictated by the 2d conformal field 
theory, should be given by llq] 



(L(O)i(f)) cx 



T.t=i\OAz,T)\ 



eiiz.r) 



El 



(8) 



where i is the temporal direction, L{x) = 
n„=i,...Ar, U^+in-i)Ls+ni is the Polyakov loop 
along i and 0^{z, t) denotes the Jacobi theta func- 



tions of argument z and modulus t = ; a 
point X — (xi, X2) in the spatial plane is mapped 
to a complex number through z = xi+ TX2 ■ Note 
that {L{0)L{x)) ~j_o , as expected at the 

critical point in the infinite volume limit. 

If the matter field behaves as an irrelevant per- 
turbation, this property should be valid even at 
/?/ > at an appropriate value of /3g(/3/). This 
has been checked accurately on large lattices at 
Nt = 6 and 20 < N^, Ny < 160. 

A typical fit is reported in Fig.(|^). Note that 
Eq.(^), being a formula derived in the context of 
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the conformal field theory, is vahd only in the 
continuum limit, so at short distance it is ex- 
pected that Polyakov loop correlators may be 
affected by lattice artifacts. Our data suggest 
that the value predicted by the continuum limit is 
reached already at distances of ~ 3 lattice spac- 
ings. Finite-size effects at criticality are rather 
strong due to scale invariance, and nontrivial. 
Therefore they are ideally suited to compare the- 
oretical predictions with MC simulations. In par- 
ticular Eq.(||) produces strong, universal shape 
effects through the modulus t, which takes into 
account the asymmetry of the lattice. A fit of the 
Polyakov correlators in an asymmetric lattice is 
reported in Fig.(||). Similar shape effects in the 
pure Z2 gauge theory at criticality where already 
observed in Ref. Il6|l. 
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Figure 8. Shape effects on the Polyakov loop cor- 
relator: same as Fig.(|^) with an asymmetric lat- 
tice of size 6 X 55 X 80. In the two figures the 
correlator is taken along the two different coordi- 
nate axes. 



Figure 7. The Polyakov loop correlator at /3g = 
0.7445, Pi = 0.17 in a lattice of size 6 x 55 x 
55 compared with the critical 2D spin spin Ising 
correlator in a square box. 

We used the goodness of the data fits to 
Eq.(||) as a criterion to locate the transition line 
Pg = fiPi) in the plane PgiPi below the SDL. 
The Kramers- Wannier transformation generates 
another critical line which is the dual /?/ = /(/?g) 
of the previous line (see Fig.(^). In the limit 
I3i -> Eq.d) yields Z(/3g, 0) cx E.,,,. ^^4z(/3g) 
where x,y,z denote periodic or anti-periodic 
boundary conditions (BC) in the 3 directions, and 
Z^y^ is the usual Ising partition function. Anti- 
periodic BC are implemented by closed surfaces 



of anti-ferromagnetic links wrapped around the 
periodic directions. 

When Pi > the sign of the links and conse- 
quently the BC become dynamical degrees of free- 
dom, however local updating algorithms on the 
critical line do not mix periodic and anti-periodic 
BC. Therefore the system behaves near the tran- 
sition line above the SDL as a pure 2D critical 
Ising model. In the Fortuin Kasteleyn (FK) ran- 
dom cluster description we can indirectly eval- 
uate the status of the BC by looking for the 
FK clusters with a linkage along the periodic di- 
rections. Transitions between periodic and anti- 
periodic BC are possible for not too large values 
of Pg- It turns out that when these kinds of tran- 
sitions become statistically relevant, the nature 
of the transition line seems modified. In partic- 
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Figure 9. The phase diagram at finite temper- 
ature. The points below the SDL are obtained 
through comparison with Eq.(H). The points 
above are obtained by comparison with the crit- 
ical partition function of 2D Ising model, as ex- 
plained in Ref.||l|]. 



ular, the agreement of the dual transition below 
the SDL worsens and the expectation value of the 
link on the transition line above the SDL is some- 
what influenced by the BC: even if the histogram 
of the distribution of the link (or the plaquette) 
variable docs not show any macroscopic double 
peak structure, we can separate this distribution 
in various sets, according to the linking properties 
of the largest FK cluster, which gives an indirect 
information on the BC of the underlying Ising 
model. The result of this separation is reported 
in Fig.(|l0|), which seems to indicate a very week 
first-order transition driven by the BC. 
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